Role of the intrinsic surface state in the decay of image states at a metal surface 
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The role of the intrinsic surface state (n = 0) in the 
decay of the first image state (n = 1) at the (111) surface 
of copper is investigated. Inelastic linewidths are evaluated 
from the knowledge of the imaginary part of the electron self- 
energy, which we compute, within the GW approximation of 
many-body theory, by going beyond a free-electron descrip- 
tion of the metal surface. Single-particle wave functions are 
obtained by solving the Schrodinger equation with a realistic 
one-dimensional model potential, and departure of the motion 
along the surface from free-electron behaviour is considered 
through the introduction of the effective mass. The decay 
of the first image state of Cu(lll) into the intrinsic surface 
state is found to result in a linewidth that represents a 40% 
of the total linewidth. The dependence of linewidths on the 
momentum of the image state parallel to the surface is also 
investigated. 



I. INTRODUCTION 

The presence of an electron in front of a solid surface 
redistributes the charge in the solid. As a consequence, 
an attractive potential is induced, which far from the sur- 
face approaches the long-range classical image potential 
Vi m (z) = —g s /£z (z being the distance from the surface, 
g s = (e s — l)/(es + 1), and e s the static bulk dielectric 
constant; for a metal, g s = 1). If the bulk band struc- 
ture projected onto the surface presents an energy gap 
near the vacuum level, an electron located in front of 
the surface cannot propagate into the solid. Therefore, 
the electron may be trapped in the vacuum well, and an 
infinite series of Rydberg-like states appears, which con- 
verges, for zero parallel momentum, towards the vacuum 
energy. These so-called image statesETB are localized in 
the vacuum region of the surface, the penetration of the 
first (n = 1) image wave function into the solid varying 
typically between 4% and 22%Q. As a result, image states 
are almost decoupled from bulk electron scattering, and 
they are much longer lived than bulk excitations: The 
lifetime of bulk electrons with energies of 4 eV above 
the Fermi level is approximately one order of magnitude 
smaller than that of the first image state. Furthermore, 
the lifetime of higher order image states [n > 1) has been 
predicted] to scale asymptotically with n 3 , which makes 
the series to be resolvable. 



During the last decade the linewidth of imag& states 
has been measured by -inverse photoemissionO'El, two- 
photon photoefflisssjoro and time-resolved two-photon 
photoemissior£3ilil. Recently, time-resolved two-photon 
photoemission has been used in combination with the 
coherent excitation of several quantum states, and the 
lifetime of the first six image states on the-(100) surface 
of copper has been accurately determined^ 

Theoretical calculations of the linewidth of image 
states were first reported in Refs. 16 and 17, within a 



many-body free-electron description of the metal sur- 
face and with the use of simplified models to approxi- 
mate both initial and final electronic states and, also, 
the screened Coulomb interaction. Later on, the decay 
of the first image state on the (111) surfaces of copper and 
nickel metals to, the n — crystal-induced surface state 
was calculated!!^, in terms of Auger transitions, with the 
use of a three band model to describe the surface band 
structure. In Ref. Q hydrogenic-like states with no pen- 
etration into the solid were used to describe the image- 
state wave functions, a simplified parametrised form was 
used for the surface-state wave functions, and screening 
effects were neglected. Self-consistent calculations of the 
linewidths of image states on copper surfaces have been 
reported recentljO, and good agreement with experi- 
mentally determined decay times has been found. In Ref. 

the linewidths of image states were computed, within 
the GW approximation of many-body theorjHI, by going 
beyond a free-electron description of the metal surface. 
Single-particle wave functions were obtained by solving 
the Schrodingerpeouation with a realistic one-dimensional 
model potentialEil, and the screened interaction was-eval- 
uated in the random-phase approximation (RPA)E3. 

In this paper we focus our attention on the role that 
the crystal-induced surface state (n = 0) plays in the 
relaxation of the first image state at the (111) surface 
of copper, which we find to represent a 40% of the to- 
tal linewidth. We present self-consistent calculations 
along the lines of Ref. [l9[ and we also consider simpli- 
fied models for both the electronic wave functions and 
the screened Coulomb interaction, showing that a de- 
tailed description of these quantities is of crucial impor- 
tance in the understanding of the origin and magnitude 
of linewidths of image states. We account for potential 
variation parallel to the surface through the introduction 
of the effective mass, and we find that the linewidth of 
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the first image state of Cu(lll) is 20% smaller than in the 
case of free-electron behavior along the surface. Finally, 
we investigate the dependence of interbandES linewidths 
on the momentum of the image state parallel to the sur- 
face, k||. Our results indicate that for image state total 
energies lying below the top of the gap, the linewidth of 
the first image state of Cu(lll) is increased with k|| 7^ 
up to a 20%. 



II. THEORY 

We assume translational invariance in the plane par- 
allel to the surface, which is taken to be normal to the 



-axis, and we evaluate the inelastic linewidth of the i 



age state </>i(z) e H' r " with energy E\ = e\ + k|/(2mi)l 

(we use atomic units throughout, i.e., e 2 = Ti — m e = 1), 
as the projection of the imaginary part of the electron 
self-energy, E(r, r';£i), over the state itself: 

r = -2 J dzdzVi(z)ImI](z,z',k||;£;i)(/.i(z / ), (1) 

where z', kii;i?i) represents the two-dimensional 
Fourier transform of X(r, r'; E*\. 

In the GW approximation! 2 !], the self-energy is ob- 
tained by just keeping the first term of the expansion 
in the screened interaction (W). Then, after replacing 
the Green function by the zero order approximation, one 
finds 
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ImX](z, z', k|| ; E\) = 

E £r E f \ E 1 

x lmW ind (z, z', q„; E 1 - E f )<j> f (z), (2) 



where the sum is extended over a complete set of fi- 
nal states </>/(z)e l ( k ii +q n)' r ii with energies £7 = e/ + 
(k|| + q|| ) 2 /(2m/). m/ is the effective mass, which ac- 
counts for the departure of the motion along the surface 
from free-electron behavior, Ep is the Fermi energy, and 
W md (z, z', q|| , E) represents the two-dimensional Fourier 
transform of the induced part of the screened interaction. 

In particular, if only transitions into the crystal- 
induced n = surface state ^>o(z)e l( - k|l+q|l ^ r|i with energy 
Eo = eo + (kii + qi|) 2 /(2m ) are considered, one finds 



ImS s (z, z', kii ; £1) 



^ 2 qn 

(2nf 



x ImVF ind (z,z',qn;£i - E )(j> {z), (3) 

and introducing this contribution to Im£(z, z', k|| ; Ei) 
into Eq. (Q) one finds: 
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dzdz'(f)l(z)^(z') 



x lmW md (z,z',q f E 1 - E Q )<f> Q {z)<t>i{z'). (4) 



Three different models have been used for the evalua- 
tion of the screened interaction, W . First, the specular- 
reflection model (SRM) of Ritchie and MarusakE^I has 
been considered. In this model, bulk electrons are as- 
sumed to be specularly reflected at the surface, the in- 
terference between the ingoing and outgoing waves being 
neglected, and the electronic charge density abruptly ter- 
minates at the surface (z = 0), which we choose to be lo- 
cated half a lattice spacing beyond the last atomic layer. 
Within this simplified modeO, also called semiclassical 
infinite barrier model (SCIBM), the screened interaction 
is obtained in terms of the wave-vector and frequency 
dependent bulk dielectric function, which we evaluate in 
the RPA. Secondly, for the vacuum contribution to the 
linewidth (z > 0, z' > 0) the surface response function 
suggested by Persson and ZarembacD (PZ) has been used. 
Finally, the screened interaction has been evaluated, as in 
Ref. [l^, by solving the RPA integral equation for the den- 
sity response function of inhomogeneous media in terms 
of the eigenfunctions of the one-electron effective Hamil- 
tonian. These eigenfunctions have been computed by 
solving the Schrodinger equation with the realistic one- 
dimensional model potential suggested in Ref. [2l|. This 
model potential uses as parameters the width and posi- 
tion of the energy gap at the f point (k|| = 0) and, also, 
the binding energies of both the n — crystal-induced 
surface state at T and the first image state. 

For the evaluation of n = and n = 1 surface states, we 
have first used simplified models for the wave functions 
inside and outside the solid. In the vacuum side of the 
surface [z > 0), n = 1 and n = states have been 
approximated by a parametrised ls-like hydrogenic wave 
function and a mere exponential, respectively, which have 
been matched to a decaying wave function in the crystal 
band-gap (z < 04-, obtained within a nearly free electron 
two-band modelE-0: 



4> n (z < 0) - e AnZ cos (Gz + 5 n ). 



(5) 



Here, n — and n = 1 correspond to the crystal-induced 
surface state and the first image state, respectively, G is 
the limit of the Brillouin zone in the direction normal to 
the surface, and 
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G 



4 gap 
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where E gap and e„ represent the energy gap and the en- 
ergy of the n surface state with respect to the midgap, 
respectively. The phase shift S n is given by 
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7r — tan 
tan -1 
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if < r] n < 1 
if — 1 < Vn < 0, 



(7) 



with T] n = 2e n /E gap . 

The image state on Cu(lll) is located right at the top 
of the gap (5i ~ 0.9 xn/2), both the hydrogenic-like wave 
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function in the vacuum [z > 0) and the decaying s-like 
wave function in the bulk (z < 0) having, therefore, nodes 
at the surface (z = 0). The n — crystal-induced surface 
state on Cu(lll) is located at the bottom of the gap (5q ~ 
0.2 x 7r/2); thus, it is described by a p-like wave function 
in the bulk. These approximate wave functions (AWF) 
are exhibited in Fig. lj together with the corresponding 
wave functions that we obtain by solving the Schrodinger 
equation with the one-dimensional potential of Ref. |2l| 
(MWF). Both wave functions, AWF and MWF, coincide 
within the bulk, but the hydrogenic-like wave function for 
the n — 1 image state appears to be less localized near the 
surface than our model wave function. The n = and 
n = 1 surface states on Cu(lll) have binding energies 
(measured with respect to the vacuum level) of 0.83 and 
5.32 eV, respectively. The n = 1 probability-density has 
a maximum at 4.3 a.u. outside the crystal edge (z = 0). 
The penetration into the bulk of n — and n — 1 surface 
states is found to be, at the f point, of 74.5% and 22.1%, 
respectively. 



III. RESULTS AND DISCUSSION 

The results of our calculations for the linewidth, r s , 
coming from the decay of the n = 1 image state into the 
n = intrinsic surface state on Cu(lll) are presented 
in Table |, with the momentum of the image electron 
parallel to the surface, k||, set equal to zero. Here, the 
linewidth has been split as follows: 
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where r vac , r so i and Fi n te r represent vacuum, bulk and 
interference contributions, respectively, as obtained by 
confining the integrals in Eq. (||) to either vacuum 
(z > 0,z' > 0), bulk (z < 0, z' < 0) or vacuum-bulk 
(z <0, z' >0) coordinates. First we show our full RPA 
calculations, in which the screened interaction is ob- 
tained on the basis of one-electron eigenfunctions com- 
puted from the realistic one-dimensional model poten- 
tial of Ref. Within these calculations both n = 
and n = 1 surface-state wave functions are also obtained 
from the model potentiaLaf Ref. |2l] (MWF), with either 
mo = 1 orpttto = 0.42E3T29. Within the specular reflec- 
tion modejEil and the model suggested by Persson and 
ZarembaEa for the screened interaction, we have used the 
MWF n = and n = 1 surface-state wave functions 
as well as the simplified models (AWF) described in the 
previous section, with too = 1. 

Total linewidths, T, as obtained from the decay of the 
n = 1 image state on Cu(lll) into any final state with 
energy Ef, E F < Ef < Ei, (see Eqs. (1) and (2)), are 
presented in Table II. Here, our full RPA calculations are 
shown, with all wave functions computed from the one- 
dimensional model potential of Ref. ^l] (MWF). Realistic 
values of the effective mass of final states have been con- 
sidered, according to the experiment or to ab initio band 



structure calculations. As for the n = surface state, we 
have used too = 0.42, as in Table I, and for bulk states 
we have chosen—to increase the effective mass from our 
computed valueEj of to/ = 0.22 at the bottom of the gap 
to to/ = 1 at the bottom of the valence band. 

Our full RPA calculations indicate that the decaying 
rate of the n = 1 image state into the n = crystal- 
induced surface state results, for too = 1, in a linewidth 
of 16 meV, while use of the more realistic effective mass 
too = 0.42 leads to a linewidth of 12 meV. With the use 
of either the free-electron mass or more realistic effective 
masses for both bulk and crystal- induced surface states, 
T s approximately represents a 40% of the total linewidth, 
T = 37meV (to/ = 1) or T = 29meV (to/ ji 1). The 
more realistic result of V = 29 meV for the total linewidth 
is in good agreement with .thr experimentally measured 
lifetime^ of 22±3 fs at 25 REE! Within the vacuum side 
of the surface the n = 1 image state couples dominantly 
to the n — surface state (this coupling approximately 
represents a 90% of the total r vac linewidth); however, 
the coupling of image states with all bulk crystal states 
occurring through the bulk penetration plays an impor- 
tant role and cannot, therefore, be neglected if one is to 
accurately describe the lifetime of image states. 

We note that simplified jcllium models for the evalua- 
tion of the screened interaction lead to unrealistic results 
for the contribution of the surface state to the linewidth 
of image statesEJ. First, we compare our full RPA calcu- 
lations (see Table I) with the results we obtain, also with 
use of our model initial and final wave functions (MWF), 
when our realistic screened interaction is replaced by that 
obtained within the specular reflection model (SRM) and 
the model of Persson and Zaremba (PZ). Bulk contribu- 
tions to the linewidth are approximately well described 
within the specular reflection model, small differences 
resulting from an approximate description, within this 
model, of the so-called begrenzung effects. As the approx- 
imate treatment of Ritchie and MarusakEa ignores the 
quantum mechanical details of the surface, this model 
fails to describe both vacuum and interference contri- 
butions to the linewidth. These quantum mechanical 
details of the surface are approximately taken into ac- 
count within the jellium model of Ref. |27|, thus result- 
ing in a better approximation for the vacuum contri- 
bution to the linewidth. Discrepancies between vacuum 
contributions obtained within this model (PZ) and our 
more realistic full RPA calculations^! appear as a result 
of the jellium model of Ref. being accurate provided 
q\\/<lF and ui/Ep << 1 (qp is the Fermi momentum, i.e., 
Ef = ql/2). 

In order to investigate the dependence of T s on the de- 
tails of both n = 1 and n = wave functions, we present 
in Table I calculations, within SRM and PZ models for 
the screened interaction, in which our realistic wave func- 
tions are replaced by the simplified models (AWF) de- 
scribed in the previous section. As the hydrogenic-like 
wave function used to describe the n = 1 image state on 
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the vacuum side of the surface presents an image-state 
charge gravity center localized further away from the sur- 
face than our more realistic model wave function, both 
vacuum and interference contributions to the linewidth 
are largely underestimated within this approximation. 
Furthermore, we note that the linewidth is highly sensi- 
tive to the details of the image-state wave functions. This 
is a consequence of the critical behavior of the imaginary 
part of the non-local self-energy coupling points near the 
surface, as we will discuss below. 

Now we focus on our full RPA calculation of the total 
linewidth of the n = 1 image state on Cu(lll) (see Table 
II), with all effective masses set equal to the free-electron 
mass. We show in Fig. ||b separate contributions to the 
linewidth, T, coming from the decay into the various / 
bulk crystal states, Tf, such that 

r = ^r / + r s . (9) 

/ 

Fig. ||a exhibits the bulk band structure projected onto 
the (111) surface of copper. The arrows indicate the 
available phase space in the decay of the n = 1 image 
state at the T point (k||=0) into the unoccupied portion 
of the n = surface state and a generic / bulk state, 
which are represented by their characteristic eo + <?|/2 

and ef + q?,/2 parabolic dispersions, respectively. r vac , 
r so i and — Tintor contributions to Tf are represented, to- 
gether with the total contribution, Tf, as a function of 

The lower edge, at the T point, of the energy gap 
projected onto the Cu(lll) surface lies below the Fermi 
level (E g < E F ) and, consequently, the decay from the 
(k|| = 0) image state occurs through finite parallel mo- 
mentum transfer. Hence, as the coupling of the image 
state with the crystal occurring through the tails of bulk 
states outside the crystal is expected to be dominated 
by vertical transitions (gy ~ 0), vacuum contributions to 
the Tf linewidth are very small, especially for those bulk 
states located at the bottom of the valence band (decay 
into these states is only allowed for large values of the 
momentum transfer; also, their vacuum penetration is 
small). Actually, the coupling of the n = 1 image state 
with all bulk states taking place at the vacuum side re- 
sults in a linewidth of only 5 meV, which approximately 
represents a 10% of the total r vac linewidthEZI. 

A realistic description of motion along the surface can 
be approximated by introduction of the effective mass, 
as described above. The effective mass of all final states 
with energies Ef, Ef < Ef < E\, is found to be smaller 
than the free-electron mass, thus both T s and T being 
about 20% smaller than in the case of free-electron be- 
havior along the surface (see Tables I and II). This is 
the result of two competing effects: First, there is the 
effect of the decrease of the available phase space, which 
is easily found to scale as yjm. Secondly, as the effective 
mass decreases the decay from the image state occurs, 
for a given energy transfer, through smaller parallel mo- 



mentum transfer, which may result in both enlarged and 
diminished screened interactions, depending on momen- 
tum and energy transfers. This is illustrated in Fig. 3, 
where the impact of the introduction of the effective mass 
on the evaluation of the various contributions (vacuum, 
bulk, and interference) to both Tf and T s (see Eq. (9)) 
is exhibited through the percentage ratio 

*r,=100^q, (10) 

as a function of the final state energy ef . If the screened 
interaction were independent of the parallel momentum 
transfer, all ratios would scale as \fm, which is repre- 
sented in Fig. 3 as R& q ,, = lOO^/mJ. Instead, as the 
parallel momentum transfer decreases the screened in- 
teraction is predominantly larger, which results in the 
ratio Rr f to be larger than RAq« > especially in the case 
of vacuum contributions to the linewidth which are ex- 
pected to be dominated by vertical transitions. When 
the decay from the image state may occur through very 
small parallel momentum transfer qu , [this is the case of 
final states that are just below the bottom of the gap and 
also the case of the n — surface state], a decrease in q\\ 
may result in a diminished screened interaction (see Fig. 
4), thus Rp f being slightly smaller than RAq« for these 
states. 

Now, we analyze the behavior of the imaginary part of 
the image-electron self-energy, which results in vacuum 
and interference contributions to the linewidth to be com- 
parable in magnitude and opposite in sign (see Table |) . 
Fig. | shows full RPA (solid line) and SRM (dashed- 
dotted line) calculations of -ImS s (z, z' , k|| = 0, Ei) (see 
Eq. (^)) for Cu(lll) with all effective masses set equal 
to the free-electron mass, together with the n = 1 image- 
state wave function, as a function of the z'-coordinate 
and for z = 7.4 a. u.. We note that the probability for 
electron-hole pair creation (the dominant channel for the 
decay of these states is provided by this process) is un- 
derestimated within the SRM. 

The coupling between electronic states is well known to 
be maximum, within the bulk, at the position of the elec- 
tron. Nevertheless, as the electron moves into the vac- 
uum the maximum of the imaginary part of the electron 
self-energy stays (see Fig. ^) near the surface (z = 0), as 
demonstrated, within a jellium model of the surface, by 
Deisz et clIl3. Hence, for any given value of z > 0, main 
vacuum and interference contributions to the linewidth 
are determined (see Eq. (|l|)) by the specific shape of 
the image-state wave function in regions A and B of Fig. 
|. As the image-state on Cu(lll) is located right at 
the top of the energy gap and the corresponding wave 
function has, therefore, a node at the surface (z ~ 0), 
an inspection of Fig. [| leads us to the conclusion that 
vacuum and interference contributions to the linewidth 
are comparable in magnitude and opposite in sign. On 
Cu(100) the image state is located close to the center 
of the gap, and the corresponding wave function has a 



4 



node at z ~ 1.3 a.u. On this surface, total vacuum and 
interference contributions to the linewidth are still op- 
posite in sign, though interference contributions coming 
from the decay into states at the bottom of the valence 
band are now positive due to their minor vacuum pen- 
etration. If the image state were located at the bottom 
of the gap, matching at the surface would occur at max- 
imum amplitude and the total interference contribution 
might be positive, as the sign of the image-state wave 
function would be the same on both sides just around 
the surface (z = 0). r— , 
Finally, we investigate the dependence of interbandE2l 
linewidths on the momentum of the image electron par- 
allel to the surface, ky. First, we use the MWF image 
state wave function evaluated at the T point and intro- 
duce, within this model, the dependence on ky of the 
quasiparticle self-energy, as indicated in Eq. (2). We 
find that for image state total energies lying below the 
top of the gap (fey ~ 0.11 a.u.), the linewidth of the first 
image state of Cu(lll) increases less than 4%. Secondly, 
we also account for the change of the z-dependent ini- 
tial wave function along the dispersion curve of the im- 
age state, by solving the Schrodinger equation for a one- 
dimensional model potential that we build following Ref. 
pl| with various values of fey: 0.06, 0.09, and 0.10 a.u.. 
The penetration into the bulk of the n = 1 image state 
with these values of feii is found to vary from 22.1% at 
fey = to 22.6%, 24.1% and 26.2% at k\\ = 0.06 a.u., 
0.09 a.u. and 0.10 a.u., respectively. Our results, as 
obtained with all effective masses set equal to the free- 
electron mass, are presented in Table III. Though the 
penetration of the image state wave function increases 
with fen , the amplitude of this wave function on the bulk 
side and near the jellium edge decreases, thus the abso- 
lute value of both bulk and interference contributions to 
the linewidth decreasing with fcii. Nevertheless, the to- 
tal overlap between image state and final wave functions 
becomes more efficient as fey increases, which results in 
larger values of the total linewidth. Also, it is interesting 
to notice that especially sensitive to the variation of the 
momentum of the image state parallel to the surface is 
the contribution to the linewidth from damping into the 
n = surface state. 



portance in the understanding of the origin of linewidths 
of image states. We have accounted for potential varia- 
tion parallel to the surface through the introduction of 
the effective mass. 

We have analyzed the origin and magnitude of the vari- 
ous contributions to the linewidth. Though the dominant 
contribution to the decay of the first image state into the 
crystal-induced surface state comes from the coupling be- 
tween image and surface states within the vacuum part 
of the surface, it appears to be approximately canceled 
out by the contribution from the interference between 
bulk and vacuum coordinates. For the vacuum contribu- 
tion to the decaying rate into the intrinsic surface state, 
we have found that it approximately represents a 90% 
of the total vacuum contribution. We also conclude that 
the coupling of image states with all bulk crystal states 
occurring through the bulk penetration plays an impor- 
tant role in the determination of lifetimes, and that this 
penetration cannot be neglected if one is to accurately 
describe the lifetime of image states. 

We have found, within our full RPA scheme, that in the 
case of ky = the decaying rate of the first image state 
on Cu(lll) into the intrinsic surface state results, with 
all effective masses set equal to the free-electron mass, in 
a linewidth of 16 meV, while use of more realistic effective 
masses leads to a linewidth of 12meV. With the use of 
either the free-electron mass or more realistic effective 
masses for both bulk and crystal- induced surface states, 
T s approximately represents a 40% of the total linewidth, 
r = 37meV (m f = 1) or T = 29meV (m f ^ 1). The 
more realistic result of V = 29 meV for the total linewidth 
is in good agreement with recent experimental results 
reported in Ref. |l3]. 

We have investigated the dependence of interband 
linewidths on the momentum of the image electron par- 
allel to the surface, showing that for image state total en- 
ergies lying below the top of the gap the linewidth of the 
first image state increases with fey up to a 20%. We con- 
clude that this increase appears mainly as a consequence 
of the change of the z-dependent initial wave function 
with fc|| , and our results indicate that the contribution to 
the linewidth from damping into the n = surface state 
is responsible for the dependence of the total linewidth 
with the momentum parallel to the surface. 



IV. SUMMARY 

We have reported calculations of the inelastic broaden- 
ing of the first image state at the (111) surface of copper, 
and we have investigated, in particular, the role that the 
intrinsic crystal-induced surface state plays in the decay 
of this image state. We have presented self-consistent 
RPA calculations, by going beyond a free-electron de- 
scription of the metal surface. We have also considered 
simplified models for both the electron wave functions 
and the screened Coulomb interaction, showing that a 
detailed description of these quantities is of crucial im- 
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TABLE I. Calculated line-width (in meV) coming from the 
decay of the n = 1 image state on Cu(lll) into the n = 
intrinsic surface state, as obtained within three different mod- 
els for the description of the surface response and two differ- 
ent models for the description of both initial and final wave 
functions (see text). RPA accounts for our full RPA realistic 
calculation of the screened interaction, SRM for the .simpli- 
fied specular reflection model of Ritchie and MarusakElj, and 
PZ for the vacuus^ side surface response suggested by Pers- 
son and Zarembao. MWF accounts for the wave functions 
obtained by solving the Schrodinger equation with the real- 
istic one-dimensional model potential of Ref. ^l|, and AWF 
for the approximate model described in the text. The effec- 
tive mass of the n = 1 image state has been set equal to the 
free-electron mass, and for the n — surface state we have 
used either mo = 1 or mo = 0.42. The momentum of the 
image electron parallel to the surface is set equal to zero. 

Surf. Res. Wave function r vac r so i ri ntcr F s 
RPA (m = 1) MWF 42 16 ^42 HT" 
RPA (m ^ 1) MWF 29 8 -25 12 



TABLE III. Calculated total linewidth (in meV) of the 
first image state on Cu(lll), computed within our full 
RPA scheme with all wave functions computed from the 
one- dimensional model potential of Ref. ^l], as a function of 
the momentum of the image electron parallel to the surface, 
fey (see text). All effective masses have been set equal to the 
free-electron mass. Contributions to the linewidth from decay 
into the n — crystal-induced surface state, F s , are displayed 
in parenthesis. 



fey 


Tvac 


r so i 


Tinter 


r 


0.0000 


47 (42) 


44 (16) 


-54 (-42) 


37 (16) 


0.0570 


48 (44) 


40 (13) 


-49 (-38) 


39 (19) 


0.0912 


50 (45) 


32 (8) 


-38 (-29) 


44 (24) 


0.1026 


50 (44) 


28 (6) 


-31 (-23) 


47 (27) 



SRM (m = 1) MWF 11 12 -17 6 

SRM (m = 1) AWF 2 15 -9 8 

PZ (m = 1) MWF 55 - - - 

PZ (m = 1) AWF 12 - - - 



TABLE II. Calculated total linewidth (in meV) coming 
from the decay of the n = 1 image state on Cu(lll) into 
any unoccupied final state with Ef < E\, as obtained within 
our full RPA scheme with all wave functions computed from 
the one-dimensional model potential of Ref. As in Ta- 
ble I, the effective mass and the momentum parallel to the 
surface of the n = 1 image state have been set equal to the 
free-electron mass and equal to zero, respectively. As for fi- 
nal (bulk and intrinsic surface) states, we have used either 
rrif = 1 or the realistic effective masses (m/ = 1) described 
in the text. Contributions to the linewidth from decay into 



bulk states, F — 


T s , are 


displayed in 


parentheses. 




Surf. Res. 


r 

vac 


r so i 


Pinter 


r 


RPA(m/ = 1) 


47(5) 


44(28) 


-54(-12) 


37(21) 


RPA(m/ / 1) 


34(5) 


32(24) 


-37(-12) 


29(17) 
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FIG. 1. Wave functions of both (a) n = 1 image and (b) 
n — intrinsic surface states on Cu(lll), as obtained within 
two different models: MWF (solid line) and AWF (dotted 
line) (see text). Full circles represent the atomic positions in 
the (111) direction. The geometrical electronic edge (z = 0) 
has been chosen to be located half an interlayer spacing be- 
yond the last atomic layer. Notice the s-like and p-like char- 
acters of the image and intrinsic-surface state wave functions, 
respectively. 
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FIG. 2. (a) Electronic surface band structure at the (111) 
surface of copper, (b) Vacuum (r vac ), bulk (F so i), and inter- 
ference (— Tintor) contributions to the linewidth (F/) coming 
from the decay into the various / bulk crystal states. The 
arrows determine the available phase space in the decay from 
the r point of the n — 1 image state into the unoccupied por- 
tion of both the n = surface state and the generic / bulk 
state. Dispersion curves of these final states are depicted. 
The energy is measured with respect to the bottom of the 
valence band. All effective masses have been set equal to the 
free-electron mass. 
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FIG. 3. Ratio Rr f of Eq. (10), as a function of the final 
state energy. The effective mass and the momentum parallel 
to the surface of the n — 1 image state have been set equal 
to the free-electron mass and equal to zero, respectively. As 
for final (bulk and intrinsic surface) states, we have used the 
realistic effective masses (m/ 7^ 1) described in the text. Here, 
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FIG. 4. Imaginary part of bulk (solid lines) and sur- 
face (dotted lines) response functions, Im [— l/e(q, u>)] and 
Im [— g s (q\\ , a>)] , which describe the screened interaction far 
from the surface into the bulk and into the vacuum, respec- 
tively. Bulk response functions have been evaluated in the 
RPA, whereas surface response functions have been evaluated] 
within the specular reflection model of the spsface (SRM)Efl 
with the RPA for the bulk dielectric functionEl 
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FIG. 5. Full RPA (solid line) and SRM (dashed-dotted 
line) calculations of — ImS s (^, z r , kii = 0,-Ei), as obtained 
from Eq. Q and as a function of the z coordinate. The 
dotted line represents the n = 1 image state wave function 
(MWF). The value of z (z = 7.4 a.u.) is indicated by an 
open circle. The sign appearing in regions A and B accounts 
for the sign of the product between the imaginary part of the 
self-energy and the image state wave function. Full circles 
represent the atomic positions in the (111) direction. The 
geometrical electronic edge [z = 0) has been chosen to be lo- 
cated half an interlayer spacing beyond the last atomic layer. 
All effective masses have been set equal to the free-electron 
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